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Abstract. We study the ageing properties of the semi-infinite kinetic spherical model 
at the critical point and in the ordered low-temperature phase, both for Dirichlet 
and Neumann boundary conditions. The surface fluctuation-dissipation ratio and the 
scaling functions of two-time surface correlation and response functions are determined 
explicitly in the dynamical scaling regime. In the low-temperature phase our results 
show that for the case of Dirichlet boundary conditions the value of the non-equilibrium 
surface exponent b± differs from the usual bulk value of systems undergoing phase 
ordering. 



1. Introduction 

Ageing phenomena encountered in systems with slow degrees of freedom are due to 
relaxation processes that depend on the thermal history of the sample. Ageing is found 
in various systems, ranging from structural glasses and spin glasses to colloids and 
polymers PU 121 E]- I n additon, ageing phenomena are also observed in ferromagnets 
quenched to or below their critical point [U EJ El E] . 

Ageing processes are best revealed through the study of two-time quantities like 
dynamical correlation and response functions. Well-known examples are given by the 
autocorrelation and autoresponse functions 



where <p(t) is the time-dependent order parameter and h is the field conjugate to <fi. The 
time t elapsed since the quench is usually called observation time and s is the waiting 
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time. In out-of-equilibrium systems one often observes in the regime t, s, t — s t micro 
{tmicro being a microscopical time scale) the following dynamical scaling behaviour [3|: 

C(t, s) = s- b f c (t/s) , R(t, s) = s-^fnit/s). (2) 

Here a and b are non-equilibrium exponents, whereas fc(y) and fuiy) are scaling 
functions which for y ^> 1 show the power-law behaviours 

fc(y) ~ y~ Xc/z , fn(y) ~ iT w * (3) 

where Ac and Xr are called autocorrelation jHUHl and autoresponse exponents [THj . and 
2; is the dynamical exponent. These simple scaling forms are observed in systems which 
are characterised by the existence of a time- dependent length scale £ ~ t 1//z . 

Specifically, for ferromagnets we distinguish between physically very different 
situations. Below T c phase ordering takes place, and the slow degrees of freedom 
are provided by the moving interfaces separating different ordered domains. We then 
commonly have 6 = 0, whereas a takes on the value 1/z or (d — 2 + T])/z, depending on 
whether the static correlations decay exponentially or whether they follow a power-law 
with an exponent d — 2 + r\ where d is the number of space dimensions. At 

the critical point general scaling arguments lead to the relation a = b = (d — 2 + r))/z 
where 7] is the well-known static critical exponent. Furthermore, the autocorrelation 
exponent Ac can be related to the initial slip exponent ^3] describing the relaxation of 
the order parameter in the short-time regime jHj. For initial states with short-range 
correlations one finds Ac = Xr, but this relation may be broken when starting from 
long-range correlated initial states [10J or when some randomness is present |15j . 

The phenomenology just briefly described has been found to be valid in many bulk 
systems. It has only been realised recently ^Hj that in semi-infinite critical systems 
similar dynamical scaling is also observed for surface quantities. At the surface we can 
define, in complete analogy with a bulk system, surface autocorrelation and autoresponse 
functions: 



c l {t,s) = (Mt)Ms)) , Ri(t,s)- ' ){o,lf)) 



ShUs) 



(t > s) (4) 

hi=0 



where <fii(t) is now the surface order parameter and hi is a field acting solely on the 
surface layer. In the regime t, s, t — s 3> t micro we expect simple scaling forms: 

C 1 (t,s)=s- b ^f Cl (t/s) (5) 
R 1 (t,s) = s- 1 -^f Rl (t/s) (6) 

and the scaling functions fc x (y) and fa (y) should display a power-law behaviour in the 
limit y — > oo: 

fcM~y- XcJz , f R Ay)~y- XRjz - (7) 
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Equations © - Q define the surface exponents ax, bx, Xd an d A^. In general 
these surface non-equilibrium exponents may take on values that differ from their bulk 
counterparts. This is similar to what is known for static critical exponents where the 
surface values differ from the bulk ones fT7\ I18j. Relations between the different non- 
equilibrium exponents can again be derived from general scaling arguments [TT)|I7]. Thus 
we obtain 



where the static exponent 7711 governs the decay of correlations parallel to the surface. 
For the surface autocorrelation exponent Aci = A^ one finds l2Ti] 



The scaling laws © - © and the different relations between the non-equilibrium critical 
exponents have been verified in ^Bj for the critical semi-infinite Ising models in two and 
three dimensions. In [7] the scaling forms in the presence of a surface have been discussed 
in more detail and calculations within the Gaussian model have been presented. 

One of the most intriguing aspects of surface critical phenomena is the presence 
of different surface universality classes for a given bulk universality class CHj- In 
this paper we study ageing phenomena in semi-infinite spherical models with Dirichlet 
boundary conditions (corresponding to the so-called ordinary transition where the bulk 
alone is critical) and with Neumann boundary conditions (corresponding to the special 
transition point where both bulk and surface are critical). Besides investigating the 
surface out-of-equilibrium dynamics at the critical point we also analyse the dynamical 
behaviour close to a surface in the ordered phase. To our knowledge this is the first 
study of surface ageing phenomena in systems where phase ordering takes place. As we 
shall see we thereby obtain for Dirichlet boundary conditions the unexpected result that 
the value of the non-equilibrium exponent bx, describing the dynamical scaling of the 
surface autocorrelation ©, is different from zero, the value encountered in bulk systems 
undergoing phase ordering 

The present work is on the one hand meant to close a gap in the study of 
kinetic spherical models [3T], which up to now have been restricted to systems with 
periodic boundary conditions. On the other hand our intention is also to extend 
towards dynamics earlier investigations of the surface criticality of the spherical model 
[2U 1221 12H1 12H 123 12B| • It has to be noted in this context that the static properties of the 
critical semi-infinite spherical model with one spherical field have been shown to differ 
from those of the O(N) model with N — > 00 [231^, even so both models are strictly 
equivalent in the bulk system [23 EHj . 

The paper is organised in the following way. In Section 2 we present the kinetic 
model in a quite general way, thus leaving open the possibility to consider different 
boundary conditions in the different space directions. In Sections 3 and 4 we discuss the 



ax = bx = (d - 2 + t)\\)/z 



(8) 



Aci = A c + 77|| - 77. 



(9) 
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dynamical scaling behaviour of surface correlation and response functions, whereas in 
Section 5 we compute the surface fluctuation-dissipation ratio. Finally, Section 6 gives 
our conclusions. 

2. The model 

2.1. General setting 

We consider a finite hypercubic system A in d dimensions containing Af — L\ x • • • x L d 
sites where L v denotes the length of the u-th edge §. To every lattice site r T = 
(ri,...,r d ) we associate a time-dependent classical spin variable S(r,t). These spin 
variables are subjected to the mean spherical constraint 

J2(S\r,t))=M (10) 

reA 

where the brackets indicate an average over the thermal noise. 

As we are interested in the rf-dimensional slap geometry we impose periodic 
boundary conditions in all but one space direction. We then have in the v-th direction 
(with v — 2, • • • , d) 

S((r 1 , ...,r u + mL u , . . .,r d ),t) = 5"((ri, ...,r„,.. .,r d ),t) for all m e Z. (11) 

In the remaining direction we either consider Dirichlet boundary conditions, 
corresponding to the ordinary transition, or Neumann boundary conditions, 
corresponding to the special transition point. For Dirichlet boundary conditions we 
impose that 

S((0, r 2 , . . . , r d ),t) = + 1, r 2 , . . . , r d ),t) = 0. (12) 

For Neumann boundary conditions, on the other hand, we have: 
S((p,r 2 ,...,r d ),t) = S((l,r 2 ,...,r d ),t) 

S({Li + 1, r 2 , . . . , r d ), t) = S{{L u r 2 , ...,r d ),t). (13) 

In the following quantities depending on the chosen boundary condition will be labeled 
by the superscript (p), (d), and (n) for periodic, Dirichlet, and Neumann boundary 
conditions, respectively. 

The Hamiltonian can be written in the following general form 

H = - l -JS T K ■ Qi T) • S A + l\ {T \t)S T A -S A -S T A -h A (14) 

where the vector S A := {S(r)\r e A} characterizes the state of the system. J > is 
the strength of the ferromagnetic nearest neighbour couplings and is chosen to be equal 
§ We set the lattice spacing equal to one. 
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to one in the following. Furthermore, h A := {h(r)\r G A} where h(r) is an external 
field acting on the spin at site r, whereas X(t) is the Lagrange multiplier ensuring the 
constraint Finally, the interaction matrix (with r = (n, • • • , r^) characterizing 
the boundary conditions in the d different directions) is given by the tensor product [21] 

d 

Qi T) = (g)(A(-) + 2E ; ,). (15) 

Here is the unit matrix of dimension L u , whereas is the discrete Laplacian in 
the z/-direction which depends on the boundary condition. 

It is convenient to parametrise the Langrange multiplier in the following way |29j : 

X^(t)=^\n)+z^(t) 

where /j a \k) is the largest eigenvalue of the interaction matrix (see Appendix A). 
In case of periodic boundary conditions in all directions we recover the usual expression 

A (p) (t) = 2d + z (p) (t). 



2.2. Langevin equation 

In order to study the out-of-equilibrium dynamical behaviour of the semi-infinite kinetic 
spherical model we prepare the system at time t = in a fully disordered infinite 
temperature equilibrium state with vanishing magnetisation. We then bring the system 
in contact with a thermal bath at a given temperature T and monitor its temporal 
evolution. Assuming purely relaxational dynamics (i.e. model A dynamics), the 
dynamics of the system is given by the following stochastic Langevin equation: 

j t S A (t) = (Q« - *M(t) - ^\n))S A (t) + h A (t) + rj A (t) (16) 

where rj A (t) := {rj r (t)\r G A} with r] r (t) being a Gaussian white noise: 

(rj r (t)) = and (7] r (t)r] r ,(t')) = 2T5 r y5(t - t'). 

We further assume that the external field h A is time- dependent. 
Equation (|TH|) can easily be solved, yielding 

S A (t) = exp (^dt'(Qi r) -z^\t') -/4 t) (k))) x [Sa(0) 

dt"(Qi T) - - fi\ K )) \ j (h A (f) + r) A (t'))] (17) 

In order to proceed further we use a Fourier-transform-like operation by defining 




t> 



o 



d 



reA \u=l 
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where the vectors \r v , k u ) are given in Appendix A. We get the original vector back 
by the inverse transformation 



fcgA \u=l 



%*)=E IK (r 1/ ,k v ))S(k,t). (19) 



The transformation (|T8|) diagonalizes the symmetric matrix Q,^ — z(t)—fJj[ (k), yielding 



--U)( T )(ft)t 



S(k,t) = -. 



S(k, 0) + jT dt' (J^it^e^^'ihik, t') + 0) 



(20) 



with 



^ (r) (t) :=exp ^2 J duz {r) [ 



u) (21) 



and 



W M(fc):=-^(fe) + ^( K ). (22) 



We now take the limit of the semi-infinite system which extends from — oo to oo in 
the directions parallel to the surface, whereas in the remaining direction the coordinate 
ri takes on only positive values. In order to stress the existence of the special direction 
we set r T = (r,x T ) with r := r\ and x T := (r 2 , ■ ■ ■ , 7"d). The corresponding vector in 
reciprocal space is then written as k T = (k, q T ) with k := ki and q T := (fc 2 , . . . , k d ). As 
a consequence of the semi-infinite volume limit sums have to be replaced by integrals in 
the following way: 

L\ IT 

• (• • •) — * f / • •) in the direction perpendicular to the surface, 
1 fci=i 71 o 

• j- Y (• • •) — > 2~ / dk(. . .) in all other directions, 

k v = l — 7T 

where the integration limits follow from Appendix A. It has to be noted that in the 
semi-infinite volume limit the largest eigenvalue of the interaction matrix (fTK|) is given 
by = 2d irrespective of the chosen boundary condition. We therefore rewrite 

(122(1 in the following way: 



u {T \k) =u(k,q) =u(k)+u(q) (23) 

with 

d 

uu(k) = 2(1- cos k) and u(q) = 2^ (1 - cos k v ) . (24) 



u=2 



where we used the explicit expressions for n^\k) (see Appendix A) 
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As usual we prepare the system at time t = in a completely uncorrelated initial 
state with vanishing magnetization. We then have at t = in reciprocal space 

(S(k, 0)S(k', 0)> = (2n) d -^5 d - 1 (q + q')C^(k, k>) (25) 

where the quantity C(k, k') is given by 

C (r) (k k>) _ f E~ i Mrk) sm(rk') r = d 

E~i cos((r - I)jfc) cos((r - \)k>) r = n W 

When we transform this expression back into direct space it just gives a decorrelated 
intial correlator, as is easily checked. 

Before deriving the scaling functions of dynamical two-point functions we have 
to pause briefly in order to discuss possible implementations of the mean spherical 
constraint in the semi-infinite geometry. Our starting point is Eq. (|10|1 which should 
hold true at all times. However, in the semi-infinite geometry translation invariance in 
the direction perpendicular to the surface is broken. One can therefore not assume a 
priori that the spins in different layers should be treated on an equal footing. In the 
past investigations of the static properties of the semi-infinite spherical model either 
considered one global spherical field [23 EH] or introduced besides the bulk spherical 
field an additional spherical field for the surface layer [2U EE] • Both cases belong to 
universality classes which differ from that of the O(N) model with iV — > oo. In 
studies of static quantities it has been proposed that the spherical constraint could 
be realised in the semi-infinite geometry by imposing in the Hamiltonian a different 
spherical field for every layer [21] , thus yielding the universality class of the semi-infinite 
O(N) model with N — > oo. For dynamical quantities, however, this prescription leads 
to the unwanted effect that the matrices appearing in the exponential in Eq. (fTTj) do not 
commute anymore, thus making an analytical treatment of the dynamical properties of 
the semi- infinite model prohibitively difficult. We therefore made the choice to retain 
only one spherical field, similar to what is done in the bulk system. This then yields 
the Equations ([27?)) - ([2~2~|) that are at the centre of our considerations. It has to be 
stressed that this implementation of the mean spherical constraint is an integral part of 
the model studied in this paper. 

3. The surface autocorrelation function 

The two-time spin-spin correlation function is defined by 

CW(r,r';t,s) = (S(r,t)S(r , ,s)) (27) 

with r T = (r, x T ) and r' T = (r, x' T ). Here (r) again characterises the chosen boundary 
condition. In the semi-infinite system correlation functions are expected to behave 
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differently close to the surface and inside the bulk. Of special interest is the surface 
autocorrelation 



C[ T \t,s) = C^((l,x),(l,x);t,s) 



(2f 



where we follow the convention to attach a subscript 1 to surface related quantities. 

From Eq. (J2Uj) we directly obtain the following expression for the correlation 
function in reciprocal space 



-u(k,q)t-uj(k' ,q')s 

C (T \k, k'; t, s) = ^(2 7 r) d - 1 -^- 1 (q + q') 

+ 2T5(k-k') [ dt'e Mk ' q)t 'g (r \t') 



C {T \k, k') 



(29) 



where we have set the external field to zero. Transforming back to real space we obtain 



C^{(r,x),(r',x);t,s) 



1 



f(r), r y. 



, t + s 



'g<T)(t)g<r)(s 
+ IT j^duf^Lr'-^-u^g^ 

where a boundary dependent function has been defined, with 



f (d) (r,r';t) 



d^q 



7T J (27r)( d - 1 ) 

f d£;sin(r- fc) sin(r' ■ k)e~ 2 ^ k)t 
Jo 

(e- 4i / (4t)) d - 1 • e- At (I r ^,(4t) - W(4t)) 



(30) 



(31) 



(32) 



and 



f in \r,r';t) = - 



d d-i q 



7T J (27r)( d - 1 ) 

p TV 

/ dk cos((r - 1/2) • k) cos((r' - 1/2) • k)e~ 2 ^ k)t 
Jo 

= (e- At I (U)) d - 1 ■ e- 4 *(W(4t) + W-i(4t)) (33) 
with the modified Bessel functions I u (u). In the long time limit t — > oo we get 

/W(r,r';t) ss 47r(87rt) _ — r • r' (34) 
/W(r,r / ;t) t «°°2(87rt)-| (35) 
where the well-known asymptotic expansion 

z/ 2 - 1/4 



Iu{u) 



V2 



ITU 



2u 



+ 0(l/u 2 ) 



(36) 



of the Bessel function has been used. 
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The function g^ T \t) is obtained as the solution of a nonlinear Volterra equation 

g<r\t) = L lim ^ ^ T \r,r;t) + 2T j\uf^(r,r;t-u)g^(u)^ (37) 

that follows directly from the implementation of the mean spherical constraint through 
the requirement 



1 L 

r lim T ^ (T) ((^0),(r,0);M) = l (38) 



r=l 

for all times t. Inserting the expressions ()34|) and (|35|) for into Eq. f)37j) we observe 
that we end up with the known bulk equation [30J 

g{t) = fit) + 2T [ duf(t - u)g{u) (39) 
Jo 

with 

/(t » = £ £ sk""** = ^ i ' i ^ d ' < 4o > 

indepedently of whether we choose Dirichlet or Neumann boundary conditions. We can 
therefore drop in the following the superscript (r) when dealing with the function g. 
The asymptotic behavior of g has been discussed in detail in [31] . We recall these results 
in Appendix B as they are needed in the following. 

Everything is now in place for a discussion of the behavior of the two-time 
correlations in the dynamical scaling regime t, s, t - s > 1. We shall in the following 
focus on the surface autocorrelation function. Other correlations (for example the 
correlation between surface and bulk spins) can be discussed along the same line. The 
surface autocorrelation function is readily obtained from Eq. (j3*Uj) by setting r = r' = 1. 

• T < Tc'- From renormalisation arguments we know that temperature is an irrelevant 
variable below the critical point and can be set to zero, and only the first term in 
()30|) contributes in the scaling limit [ST]. For both boundary conditions we obtain 
a dynamical scaling behavior, as 

CW(t,s)=2tM? q s-i(t) (t + lj , (41) 

- \ — - 

^)(t, S ) = 2l+ 1 M e 2 g 0j 4 + lj 2 . (42) 

where M^ q = 1 — T /Tc |H2]- The non-equilibrium exponents can then be read off 
directly (recall z = 2): 

b[ d) = l, &i n) = 0, Ag> = ^ + 2, \%= d -. (43) 
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A few comments are now in order. When comparing the values (J43|) with those 
obtained for the bulk system see Tabled we observe that for Neumann 

boundary conditions the non-equilibrium exponents take on exactly the same values 
than in the bulk. This is not only observed for quenches to temperatures below the 
critical point but also for quenches to the critical point, see below, and also holds 
for the exponents related to the response function, as discussed in the next Section. 
Interestingly, the agreement between Neumann and periodic boundary conditions 
(i.e. bulk systems) also extends to the scaling functions of two-time quantities 
which are identical, up to an unimportant numerical factor. | 
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Table 1. Ageing exponents and fluctuation-dissipation ratio for periodic (bulk 
system), Dirichlet and Neumann boundary conditions. 

The situation is different when considering open (i.e. Dirichlet) boundary 
conditions. The values of the non-equilibrium exponents are distinct from the values 
of their bulk counterparts and the scaling functions are different, too. The larger 
value of the surface autocorrelation exponent A[? thereby reflects the increased 

|| This behavior may be compared to that of the critical semi- infinite Ising model at the special 
transition point where the values of the non-equilibrium surface and bulk exponents and the surface 
and bulk scaling functions are found to disagree Similarly, the surface autocorrelation exponent 
Ad for the O(N) model in the limit N — ► oo differs at the special transition point from the bulk 
autocorrelation exponent Ac Differences at the special transition point between the semi- infinite 
spherical and the Ising and O(N) (N — ► oo) models already show up in the static properties as the 
surface and bulk order parameter exponents are identical for the spherical model but differ for the other 
two models [TH) . 
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disorder at the surface due to the absence of neighbouring spins. Remarkably, 
the exponent is found to be different from zero, the value usually encountered 
when studying phase ordering in a bulk system. It is an open question whether 
this surprising observation is unique to the special situation of the spherical model 
or whether it is encountered in other systems as for example the semi-infinite Ising 
model quenched to temperatures below T c . 

The two types of dynamical scaling behaviour encountered in the semi-infinite 
spherical model quenched below the critical point are illustrated in Figure ^ in 
three dimensions. The analytical curves are given by the expressions (|41|) and (|42jl 
whereas the symbols, corresponding to different waiting times, are obtained from 
the direct numerical evaluation of Equation (|30|). It is obvious that we are well 
within the dynamical scaling regime even for the smallest waiting time considered. 
Furthermore this confirms a posteriori that it was justified to drop the second term 
in flU). 




02468 02468 
t/s t/s 

Figure 1. Scaling plots of the autocorrelation function for the case T < Tc in three 
dimensions: (a) for Dirichlet boundary conditions with b± = 1, (b) for Neumann 
boundary conditions with = 0. The inset in (a) shows that no data collapse is 
observed for the unsealed data. 

• T = Tc and 2 < d < 4: The behavior of C[ (i, s) in the regime t, s, t—s ^> 1 follows 
from the insertion of the asymptotic expressions for g and for f( T > into Eq. 1)30)1 . 
One may wonder whether the use of the asymptotics for g in the integrand of Eq. 
(J30j) does not cause problems at the lower integration bound. As has been shown 
in [34J there exists for the spherical model a time scale tp ~ with < £ < 1 such 
that for times larger than tp one is well within the ageing regime. Replacing g by 
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its asymptotic value g age for u > tp, the integral of Eq. (J3UJ) can be written in the 
following way: 

f S duf T \l,l;(t + s)/2-u)g(u) 
Jo 

= I" dug(u)f T \l 7 l;(t + s)/2-u) + [ dug age (u)f^(l, 1; (t + s)/2 
Jo Jt P 

= W {T) (t,s,t P ) + s [ dvg age {sv)f {T \l,l-{t + s)/2-sv) 
Jo 

where in the last line we have assumed s to be large. An upper bound for the first 
term W^(t, s, t P ) = f* P dug(u)f^{l, 1; (t + s)/2 - u) is given by 

\W p T \t,s,t P )\ <t P max \g(r)f^(l,l;(t + s)/2-r)\ 



re[o,t P ] 



S>1 



re[o,t P ] 



tp ■ max \g(r)\s ? 47r -(- + 1) 



1 ,t 



tp 



2 K s 



where £^ = | + 1 and £ (n ) = |. As t P ~ s^, this upper bound disappears in the 



scaling limit faster than the second contribution in Eq 
W^ T \t, s, tp) in the following and end up with the expression 



(45) 
i the 

We therefore drop 



CJ T) (M) 



\/9age(t)g 



/ (r) 1,1; 



t + s 



+ 2T / duf 



e(r) 



i,i; 



t + s 



u 



g age (u)d 



a 



(46) 



in the asymptotic regime. It turns out that the thermal part of the expression (J46j) 
is the leading one at the critical points in any dimensions, so that the first term can 
be disregarded. For 2 < d < 4 we then obtain again dynamical scaling behavior, as 
shown in Figure as: 



C (d) (t,, 



4T c (4tt) 



-(f) 



d-2 



2 ' -± 
S 2 



4 ft 



C {n \t, 



d \s 
8(4tt) 



t 

- + 1 

s 



t 

s 

-2 



1-'- 



-1 



T, 



£.,-(1-1) [ t 
s 



(n) 



(d-2) 

yielding the non-equilibrium critical exponents 
bl "2' 

The values b^ agree with the expression 

2(,i?! T) - P) 



s 

i-i' 



(47) 



s 



1-1 



- + 1 (4 



d 



Ag = ~d-2. (49) 



= b + 



vz 



(50) 



expected from general scaling considerations 7J. Eq. (j50p is readily verified in three 
dimensions where b — ^ = |, whereas j3 and [3[ T \ which are the usual static 
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equilibrium critical exponents describing the vanishing of the bulk and surface 
magnetisations on approach to the critical point, take on the values /3 — \ and 
Pi = | ? Pi^ = \- Finally, the static exponent v = 1 in three dimensions. 
As a final remark, let us note that the tendency of the surface correlations to decay 
faster in time than the bulk correlations is mirrored at the ordinary transition by 
the larger value of the non-equilibrium autocorrelation exponent Ag? . 
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Figure 2. Scaling plots of the autocorrelation function for the case T = Tc in three 
dimensions: (a) at the ordinary transition (Dirichlet boundary conditions), (b) at the 
special transition point (Neumann boundary conditions). 



T — Tn and d > 4: In this case we obtain 



C {d) (t, 



C {n \t, 



2(4vr) 



d 

4(4tt) 



s 



t 

- + 1 

s 



s 



s 



(51) 



(52) 



d- 2 

Again dynamical scaling is found where the non-equilibrium critical exponents now 
take on the values 

d d 
2' 



id) 



- 1, 



(53) 



Also here (J5U|) is verified with the values 6 = ~ — 1, = 1, /3j = |, z/ = | and 

/3 



i 



2 ' 

T > Tq: We then have 



C{ T) (M) = e-( t+ ^ 27 --) 



/ (T) (M, 



t + s 



Semi-infinite kinetic spherical model 14 



+ zr jr-d»/M(i,i : *±i-«) 



U I g a ge{u) 



(54) 



This amounts to an exponential decrease and the leading expression rapidly 
develops a dependence only on t — s for s — > oo and t — s fixed [HI] . 

Inspection of Table ^ reveals the oddity that the values of a\ and b\ for Dirichlet 
and Neumann boundary conditions always exactly differ by one. This may again be 
compared with the critical semi-infinite Ising model. We there have a\ — bi — 2f3\/{vz) 
[inilZj with the values /3i = 0.80 at the ordinary transition and (3% ~ 0.23 at the special 
transition point [TH], yielding different values for a x and b\ in that case, too. However, 
the exact value one for the difference seems to be a property of the spherical model. 



4. The response function 

The response function is defined as usual: 

8(S(r,t)) 



R(r,r';t,s) := 



6h(r', s) 



(55) 



h=0 



where h(r', s) is a magnetic field acting at time s on the spin located at lattice site r' . 
Starting from expression (j2U|) and assuming spatial translation invariance parallel to the 
surface, we find 



R((k,q),(k',qy,t,s) = ^ 9 j^-e-^^5{k-k'). (56) 
The chosen boundary condition enters when transforming back to real space, yielding 



R( d \(r,x),(r>,y);t,s) = J ^ e ~ 2 ^ d (I B _ y (2(t - s)))^ 1 



x (I r ,_ r (2(t-s))-I r/+r (2(t-s))) (57) 



R^((r,x),(r',y);t,s) = ^J^ e -^ d (/ (B _ 1/ (2(t - s))Y 

x (J r /_ P (2(t-s))+/ r / +P _i(2(t-s))). (58) 

We refrain from giving a full discussion of the response function, but focus instead on 
the surface autoresponse function Ri(t, s) = R^((l, x), (1, x);t, s) that describes the 
answer of the surface at a certain position to a perturbation at the same site at an 
earlier time. With the large t behaviour (|3*B|) of the Bessel function we then find 



R^it^s) =4vr 



J^iMt - s))^ (59) 



RP(t, s) = 2j||(47r(t -*))-!. (60) 
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It is evident from these expressions that the ageing behaviour of the surface autoresponse 
is again determined by the asymptotics g age of the function g given in the Appendix B. 
As for the autocorrelation function we have to distinguish four different cases: 

• T < Tc'- Here we have 

R[ d \t,s) = (An)-il^\ s -(l+i)(I-lj (61) 

- \ — - 

4 n \t, s) = 2(4tt)-^) (t\ 1 s-l ft - 1 j \ (62) 

For both choices of boundary conditions we observe dynamical scaling (see Figure 
BJ) with the exponents 

«. M = 5, <4 W = |-1, Ag = | + a, Ag = |. (63) 




Figure 3. Scaling plots of the autoresponse function for the case T < Tc in three 
dimensions: (a) for Dirichlet boundary conditions with = |, (b) for Neumann 
boundary conditions with d[ n ^ = h. 



T = Tc and 2 < d < 4: We find again dynamical scaling, as 



>(»)' * • - ; 



t\-i+ 1 d ft \-(f+ 1 ) 



RW(t, s) = (4tt)-2 ( II s -lf+iJ 11 - ij (64) 

d i i d 



i?f ; (t,s) = 2(4tt)~5 ( I J s -ifl_i) (65) 
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Figure 4. Scaling plots of the autoresponse function for the case T = Tq in five 
dimensions: a) at the ordinary transition (Dirichlct boundary conditions), (b) at the 
special transition point (Neumann boundary conditions). 



with the following nonequilibrium critical exponents: 

Ad) _ d (n) _ d {d) _ 3 (n) _ 3 

a i - 2' 1 ~ 2 ' Rl ~ 2 Rl ~ 2 
T = Tq and d > 4: This case yields the expressions 



R[ d \t,s) = (4tt)- ( - - 1 ) (67) 



s 

R^\t, s) = 2(47r)-ls-i ^ ~ l) 2 • ( 68 ) 
Again dynamical scaling is found, the exponents now taking on the values 

J?> = i a?= d --l, Ag = d + 2, Ag = A (69) 

This dynamical scaling behaviour is illustrated in Figure 0] in five dimensions. 

T > T g : Due to the exponential behaviour of g a ge(t), R(t,s) disappears 
exponentially: 

R[ d \t,s) = 47re- ( ^ s)/re «(4vr(t - s))~^ +1) (70) 

R^it, s) = 2e~ (i - s)/Te «(47r(t - s))~5 (71) 

and only a dependence on the time difference t — s is observed. 
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Comparing the values of non-equilibrium exponents derived from the autoresponse with 
those obtained from the autocorrelation reveals that we have at criticality the identities 
off = bi and A^ = A^ in agreement with the general scaling arguments given in the 
introduction. 



It has been argued [35] that in out-of-equilibrium dynamics the scaling functions of 
response functions can be derived from the space-time symmetries of the corresponding 
noise-less Langevin equation. This has been checked for a series of exactly solvable 
systems (see [HUE] and references therein), among them the spherical model in the bulk 
system ^f. In [HI] local scale invariance was used to derive the following prediction for 
the surface autoresponse function (with t > s): 

R L l s \t,s)=rJ-\ 2 \t-s)^ (72) 



s J 

where (\ and (2 are two exponents left undetermined by the theory, and tq is a non- 
universal normalization constant. The values of £i and (2 are fixed by comparing (J72j) 
with the expected scaling behaviour (see Equations (0) and (J7J)) yielding the result 

R^ I (t,s)=r s-^ ) r-) . (73) 



This prediction is in complete agreement with our exact results for T < T c both for 
Dirichlet and Neumann boundary conditions. 

Finally, let us mention that it was shown in [3H] that the dependence of the full 
response function ()58|) on r and r' can be deduced from the space-time symmetries in 
the case of Dirichlet boundary conditions. 



5. The fluctuation-dissipation ratio 



The fluctuation-dissipation ratio 



TR(t,s) 

{,) - d s C(t,s) U4j 



has been discussed extensively in recent years j39j as a possible way for attributing an 
effective temperature to an out-of-equilibrium system. Of importance in the following 
is the fact that a characteristic behaviour is expected for different physical situations. 
Thus in systems undergoing phase ordering one usually expects X(t, s) to approach 0, 

% Note, however, that in critical systems with z ^ 2, as for example in the critical two- and three- 
dimensional Ising models, corrections to the prediction of local scale invariance (LSI) have been shown 
to exist . But even in that case the scaling function coming from the theory of local scale invariance 
|H7j yields an excellent description of the numerically determined autoresponse function. 
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whereas a different behaviour is observed for non-equilibrium critical systems. Indeed 
it has been realized [HH U\ that in the latter case the non- vanishing limit value 

X°° : = lim \im X(t,s) (75) 

is a universal quantity whose value characterises the given dynamical universality class. 

The concept of fluctuation-dissipation ratio, initially introduced for bulk systems, 
has been generalized in [TH] to systems with surfaces. The surface fluctuation-dissipation 
ratio is thereby defined by 

Xl(M) ' = ^M) (76) 

The asymptotic value X^° := lim^oo lim^oo Xi(t, s) is in fact a ratio of two amplitudes 
and its value at the bulk critical point should be characteristic for a given surface 
universality class |7j. Recently, Xi(t, s) has been determined numerically for the critical 
two- and three-dimensional semi-infinite Ising models JH], with asymptotic values X^° 
differing from the values X°° obtained in the corresponding bulk systems. X^° has also 
been computed within the Gaussian model |7j yielding the value |. 

Having already determined the surface autocorrelation and autoresponse functions 
in the previous sections, we are in the position to compute Xi(t,s) also for the semi- 
infinite spherical model. In doing so one immediately realizes that for Neumann 
boundary conditions x[ n '(t, s) is identical to the bulk quantity X(t, s) obtained in |3Tj . 
This follows from the fact that C\ and R\ for Neumann and periodic boundary conditions 
only differ by a numerical constant which drops out when the ratio is formed. For this 
reason we give here only the results obtained for Dirichlet boundary conditions: 

• T <T C : 



d 



X«>(t s) - lMZ g -f + i j±i (77) 

Xl Mi q (A +d y--d (77) 

yielding the limit value X^° = as expected for ferromagnetic systems quenched 
below their critical point. 

T = T c and 2 < d < 4: In this straightforward but somewhat tedious 

calculation yields 

X W (t8) = d{d-2){^lf 

(d 2 -16) + (32 + d(3rf-8))f + (rf-4)(4 + 3rf)(^) 2 + rf 2 (^)3 1 ] 

with the limit value Xf = 

T = Tc and d > 4: Here we obtain the same expressions as for Neumann and 
periodic boundary conditions: 

__(*^ 1 ,. 
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with X 



1 



oo 



1 

2' 



In all cases the limit value turns out to be independent of the boundary condition, which 
is the reason why we have dropped the superscript (d). We therefore conclude that in 
the semi-infinite spherical model X%° equals X°° independently of the chosen boundary 
condition. 

6. Conclusion 

In this paper we have extended the study of out-of-equilibrium dynamical properties 
of the kinetic spherical model to the semi-infinite geometry with both Dirichlet and 
Neumann boundary conditions. The exact computation of two-time surface quantities 
(like the autocorrelation and the autoresponse functions) reveal that dynamical scaling 
is also observed close to a surface for quenches to temperatures below or equal to the 
critical temperature. Whereas for Neumann boundary conditions we find that the values 
of the non-equilibrium exponents and the scaling functions (up to a numerical factor) 
are identical to the bulk ones, the situation for Dirichlet boundary conditions is more 
interesting. Indeed at the critical point we find that out-of-equilibrium dynamics is 
governed by universal exponents whose values differ from those of the corresponding 
bulk exponents. The values of these surface exponents are in complete agreement with 
predictions coming from general scaling considerations Similarly, surface scaling 

functions are also found to differ from the bulk scaling functions. Interestingly, we find 
that in the ordered low-temperature phase the autocorrelation function scales in the 
ageing regime as 



As this paper is the first study of ageing phenomena close to a surface in an ordered 
phase, it is an open und interesting question whether this is only a special feature of 
the kinetic spherical model or whether this is a general property of semi-infinite systems 
undergoing phase ordering. We intend to come back to this problem in the near future. 

Appendix A. Eigenvalues and eigenvectors of the interaction matrix 

The eigenvalues (fc) of the interaction matrix depend on the boundary 

conditions and are given by j2Hj 

d 




(80) 



in strong contrast to the well-known bulk behaviour 0] 






(Al) 
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(t) 

where the functions (p L are 

( -Pl\^) = ~~r~ (periodic) 
<pf(k) = (Dirichlet) 

{n) (u \ ^ik-l) 

L PL\' t ) = t (JNeumannJ 

LJ 

The (orthonormal) eigenvectors are given by the expressions 



u A 



with 



uf(r,k)—L 2 exp(— inp^ik)) 



u 



r,k) = [2/(L + l)]ssin(r^(Jfc)) 



U " ){r,k) ^ (2/L)lcos((r-|Vi 0) (A;)) k = 2,...,L 



Appendix B. Asymptotic behaviour of g 

In this Appendix we recall the asymptotic behaviour [31] of the function 

g{t) = fit) + 2T [ duf(t - u)g{u) 
Jo 



with 



m ■= (2^1 y^ e ~ MM< = ^%m d - 

Defining first the constants 

d d-i q «r dk 1 



{2n) d - 1 J_ n 2tt [2u{k,q)) k 



the critical temperature can be expressed as 

1 



T, 



c 



2A 1 

The asymptotic behaviour of g for long times, denoted g age , is then given by 
• T>T C : 

9age(t) e 



t/Teq 



where the exact expression for the time scale r eq can be found in [3Tj . 
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T <T C : 



0*e(*) ~ With K = 1 -tT> ( B6 ) 

T = T c and 2 < d < 4: 

fl*e(f) - (d - 2)(87r)^ 1 sin ( ) (B7) 



T = T c and d > 4: 
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